Density functional theory in one-dimension for contact-interacting fermions 
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A density functional theory is developed for fermions in one dimension, interacting via a delta- 
function. Such systems provide a natural testing ground for questions of principle, as the local 
density approximation should work well for short-ranged interactions. The exact-exchange contri- 
bution to the total energy is a local functional of the density. A local density approximation for 
correlation is obtained using perturbation theory and Bethe-Ansatz results for the one-dimensional 
contact-interacting uniform Fermi gas. The ground-state energies are calculated for two finite sys- 
tems, the analogs of Helium and of Hooke's atom. The local approximation is shown to be excellent, 
as expected. 

PACS numbers: 31.15.Ew, 71.15.Mb, 71.10.-w, 



I. INTRODUCTION 

Density functional theory (DFT) is a rigorous re- 
interpretation of the quantum many-body problem in 
which the basic object uniquely characterizing a system 
is the density, n{x), rather than the many-body wave- 
function. This view is particularly suited to practical 
calculations, and DFT has been applied successfully to 
solids and molecules for quite some time [1] . Researchers 
typically focus on Coulomb-interacting fermions in three 
dimensions, but the Hohenberg-Kohn theorem, or the 
one-to-one correspondence between potentials and den- 
sities [2] upon which DFT is based, holds for any inter- 
action and in any spatial dimension. We consider the 
contact or delta- function interaction between fermions in 
one spatial dimension, 



= (Xi -^J ) — XS(^Xi '^j^: 



(1) 



where Xi and Xj represent the spatial coordinates 
of the fermions, 6{x) is the Dirac delta function, and A 
is the interaction strength. The fermions have two spin 
states, up and down. The delta-function potential is a 
one-dimensional analog to the Coulomb one as it scales 
in a similar fashion: Vee{ax) — Vee{x)/ a, and its solutions 
satisfy the energetically important particle-particle cusp 
condition [5, 6]. However, it differs in that it is short- 
ranged. There is no simple equivalence between A and 
e^, the Coulomb- interaction strength, although A can be 
related to a scattering length [4] . 

One-dimensional model interactions are important for 
several reasons. Perhaps most obviously, they are useful 
in mathematical and statistical physics [7-9] to illustrate 
problems and concepts from three-dimensional physics 
that are sometimes hard to conceptualize due to the num- 
ber of degrees of freedom. However, our primary moti- 
vation is to use this one-dimensional model to under- 
stand and improve density functional theory. Many of 
the known formal properties of the exchange-correlation 



functional are true in this case. These properties include 
behavior under uniform coordinate scaling[10], the virial 
theorem, and inequalities due to the variational prin- 
ciple. That properties of the theory still hold should 
prove extremely useful in exploring time-dependent den- 
sity functional theory[ll, 12], where formal properties are 
still being explored. Because the interaction is not the 
Coulomb one, the unknown exchange-correlation func- 
tional will, of course, differ. The local density approxi- 
mation should be extremely accurate in this case because 
of the short-ranged nature of the interaction. Thus in our 
case, failures of the local density approximation (LDA) 
can be ascribed to non-locality that is independent of the 
long-ranged nature of the Coulomb repulsion. Our LDA 
could be used to study the one-dimensional analog of 
stretched H2 to identify whether the proper description 
of dissociation into individual atoms depends on the long- 
ranged Coulomb interaction or is due to symmetry con- 
siderations alone. Another interesting system on which 
to use one-dimensional DFT is the one-dimensional solid. 
This delta-function interaction has already been used to 
study problems in DFT [13] but without the inclusion of 
any correlation effects, which are known to be important 
in one-dimensional systems. Using DFT to study alter- 
nate interactions is not new; for example, Capelle and 
coworkers have used a similar approach on the Hubbard 
model[14, 15]. 

It has been suggested that the delta- function model 
should give a good representation of the physics of one di- 
mensional fermions in certain experimental contexts [16- 
20] . Since one-dimensional systems are analytically, or at 
least computationally, manageable, the exact results are 
useful to examine situations when standard techniques 
fail. For example, the one-dimensional analogue of He- 
lium can be examined in detail near the critical point of 
ionization when the nuclear attraction and interaction re- 
pulsion are comparable. This sort of analysis is demand- 
ing for real systems, and finite-size scaling and infinite 
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dimensional approaches are necessary [21-24]. Carefully 
understanding how systems ionize and how electronic 
structure methods reproduce this critical phenomenon 
are useful for many chemical problems. We will examine 
this limit in detail in future work and present only the 
most basic results here. 

Throughout, we assume that our one-dimensional 
fermions have the same mass as electrons, and we use 
atomic units {e^ = U = rUe = 1) so that all energies are 
in Hartrees and all lengths in Bohr radii. 



II. EXACT-EXCHANGE FUNCTIONAL 
DENSITY FUNCTIONAL THEORY 



In this section, we see how the contact interaction af- 
fects the total energy to first order in A. First order inter- 
action theory is traditionally called the Hartree-Fock ap- 
proximation, but here, the first order interaction energy 
depends explicitly on the density so that, for this par- 
ticular interaction, Hartree-Fock is equivalent to exact- 
exchange DFT. Consequently, exchange is treated ex- 
actly within the local density approximation for this in- 
teraction. 

According to the Hohenberg-Kohn theorem [2], the 
ground-state total energy is a functional of the particle 
density: 

E[n] = Ts[n] + U^n] + E^a[n\ + j dx v^^t{x)n{x) (2) 

in a one-dimensional space where J7h['t^] is the exactly 
known Hartree or classical density-density interaction 
contribution, Woxt(a^) is the given external potential, Ts[n\ 
is the exactly known kinetic energy of non-interacting 
fermions at a given density, and £^xc["^] is the un- 
known exchange-correlation energy. The density is found 
by studying the Kohn-Sham (KS) system, the non- 
interacting counterpart to the physical system [25]. The 
Kohn-Sham equation is 

(^-^'^'^ +vs,a{[ni,ni];x)j (t>i,a{x) = ei,a(l>i,a{x) (3) 

where (pi^a{x) is the i*^ KS orbital for spin-type, a, ei^a 
is the KS eigenvalue, vs,a{x) is the KS potential for spin- 
type, a, and and rij are the densities of up- and down- 
spins. The Kohn-Sham potential is a functional deriva- 
tive of the energy functionals, 

-sAx)=vMx) + ^^+^, (4) 

where na{x) is the cr-spin density. The spin-density is 
obtained from the occupied orbitals. 



na{x) = ^ \<l>i,a{x)\''' 

i.OCC. 



(5) 



Because of the anti-symmetry of the wave-function un- 
der particle interchange, fermions with like spins do not 



experience the contact interaction. Only opposite spins 
interact directly. 

The Hartree contribution depends only on the total 
particle density and is independent of how up and down 
fermions are distributed: 

UHin] ~ ^ J '^^ ''T'{x)S{x — x')n{x') 

= \jdxn\x). (6) 

There is over-counting here because like spins do not in- 
teract, and the exchange term must cancel these spurious 
like-spin interactions. The exact-exchange term is 



Ex 



[m,ni]=-'^ j dx {n]{x) + nl{x)) . 



(7) 



Both terms can be derived using the standard rules of 
many-body perturbation theory [26] (see appendix A). 
To simplify the notation, let 



ax) = 



n|(a;) - ni{x) 



Then, the exact-exchange functional can be written, 
ExM=-\ j dxn{xf{l + axf)/2, 



(8) 



(9) 



Note that for a one fermion system, we have £^jc[n] = 
— C//f[n] and contact-interacting exact-exchange is self- 
interaction free. 



III. LOCAL DENSITY CORRELATION 
FUNCTIONAL FROM DELTIUM: THE 
ONE-DIMENSIONAL FERMION GAS 

In order to obtain a local-density correlation func- 
tional, we review the one-dimensional unpolarized Fermi 
gas, which we call deltium. This Fermi gas plays the role 
of the uniform electron gas in Coulomb-interacting DFT. 
While the Coulomb-interacting Fermi-gas is a Fermi liq- 
uid, the one- dimensional delta-function interacting ana- 
log is a Luttinger liquid [3] . The Hamiltonian is 



1 ^ d2 



(10) 



The solution must be anti-symmetric under particle in- 
terchange and satisfy periodic boundary conditions on a 
ring of circumference, L. This system has been examined 
previously [4, 27, 28]. Because the wave- function is anti- 
symmetric under particle interchange, the fully polarized 
gas is not affected by the interaction. We consider the 
correlation in detail for only the fully unpolarized gas. 

The energy per particle of the non-interacting uniform 
gas is purely kinetic: 



(11) 
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When interactions are present, the total energy per par- 
ticle is 



e(n) = t{n) + eH{n) + e^{n) + ec(n) 



(12) 



where eH(n) = An/2 is the Hartree energy per particle, 
ex(n) = —An/4 is the exchange energy per particle, and 
Cc (n) is the correlation energy per particle. It is useful to 
define the following two terms. Kinetic-like means that 
the energy per particle is proportional to "n? like the non- 
interacting kinetic energy. Hartree-like means that the 
energy per particle is proportional to n like the Hartree 
energy. 

The ground-state energy per particle for deltium, Eq. 
(10), can be found via Bethe Ansatz methods [29-33]; 

whereby, the uniform unpolarizcd Fermi gas problem can 
be recast as a set of integral equations [34] : 



r{y) = — + - I dA- 



2tt tt 



A2 + 4(2/ - A)2 



(13) 



and 



^(A) = 2T / sech (7r(y - A)/A) r(y) (14) 

where r is the number of occupied states per wave- vector 
label, y; and a is the number of occupied down-spin 
states per a different wave-vector label, A. In the high- 
density limit, kmax = 7rn/2, and in the low-density limit, 
kmax — Trn. Equations (13) and (14) must be solved 
self-consistently for t and ct at a chosen value of kmax to 
obtain the ground-state energy. In order to do this, the 
integrals are transformed to the interval, [-1,1], and in- 
tegrated using six-point quadrature rules with 400 mesh 
points. The density is 



J-i 



dy T{y), 



and the total energy per particle is 



~ J-K 



dy T{y). 



(15) 



(16) 



The correlation energy per particle for a wide range of 
densities can, in the spirit of three-dimensional DFT, 
be parameterized for practical calculations. We consider 
both the high and low-density limits analytically and nu- 
merically. Since we are concerned with parameterizing 
correlation energy, we subtract out the known kinetic, 
Hartree, and exc;liange contributions. 

The low-density limit is the large A limit. In this limit, 
the opposite-spin Fermions repel each other so strongly 
that the interaction mimics Fermi anti-symmetrization. 
Thus, the interaction energy per particle is kinetic-like. 
This means that the first term for the correlation energy 
must cancel the Hartree energy, and the next term in the 



correlation energy must be kinetic-like. The correlation 
energy per particle is 



e^{n) = -'-^n + bin'^ -^n^ + 



A 



(17) 



with hi = TT^/S and 62 = 4.560971. The first terms in Eq. 
(17) exactly cancels the Hartree and exchange energies, 
and the next, b\, provides the kinetic- like contribution 
to the energy. We determine the coefficient, 62, from 
numerical analysis of the Bethe Ansatz resiilts at A = 1. 
Specifically, we subtract the known contributions from 
e(n), divide by n^, and plot the result as a function of n. 
The extrapolated intercept at n = is 62- 

In the high-density limit, the interaction is 
perturbation-like and the correlation energy per particle 
approaches a constant: 



A'' 

ec(n) = -ciA^ + C2>?/n + 0[^ 



(18) 



with ci = 1/24 and C2 = 0.006151. The first term, ci is 
found in appendix A using perturbation theory. We de- 
termine C2 from numerical analysis of the Bethe Ansatz 
results. We find the expansion coefficients by subtract- 
ing the known contributions from e(n) and plotting the 
remainder as a function of 1/n. The coefficients are the 
extrapolated intercept and slope as a function of 1/n. 
The correlation energy per particle of the uniform gas 
approaches a finite value as n ^ 00 in contrast to jel- 
lium, because the contact interaction is short-ranged. 

A [2,2] Fade parameterization of the correlation energy 
per particle is 



in) 



bn ■ 



v? + dn + e 



(19) 



with a = -0.0416667, b = 0.004475, c = 0,d = 0.254998, 
and e = 0.017900. Note that c is zero because the correla- 
tion energy per particle vanishes in the low-density limit. 
This approximation gives the first and second terms of 
both the high and low-density limits correctly. The pa- 
rameterization of Ec, the correlation energy per particle, 
has a maximiim error of 0.7% and is highly accurate for 
the important high density regions. As shown in Figures 
1 and 2, the parameterization is almost indistinguishable 
from the exact numerical result. 

For the fully polarized case, the interaction does not 
contribute and 



i;^^^[n] = 0. 



(20) 



We can combine these results and construct a local- 
density correlation energy functional: 

Eb''>] = jdx n{x) eTc'^-{n{x))mx)) (21) 
an{x)^ + bn{x) 



-hi 



n{x)'^ + dn{x) + e 



fiCix)). 



(22) 
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FIG. 2: The same as Fig 1 but for the low-density hmit. 
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TABLE I: Total ground-state energy for Diracium and er- 
rors (in milliHartrees) of various approximations with A = 1. 
The exact results are from a numerical solution of the prob- 
lem as outlined by Ref. [36], the second-order perturbation 
values axe also given in the same reference. The EXX and 
the LDA results are from a self-consistent solution using the 
exact-exchange functional, Eq. (7), and the LDA functional 
4 of Eq. (22) respectively. Zcrit = 0.377115 is the critical value 
at which the system ionizes [36]. 



FIG. 1: Correlation energy per particle for Deltium, the onc- 
dimensional uniform contact-interacting (A = 1) Fermi gas, 
in the high-density limit. The solid line is the exact result 
calculated from the solutions of the Bethe-Ansatz integral 
equations. The long-dashed line is the simple LDA param- 
eterization given by Eq. (22). 
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FIG. 3: Correlation energy for Diracium, the one-dimensional 
analog of Helium (A = 1). The solid line is the exact result, 
the short-dashed line is the LDA result, and the long-dashed 
line is the second order perturbation theory result. Note that 
we extract the correlation energy approximately by subtract- 
ing the self-consistent exact-exchange total energy from the 
exact total energy. The exact-exchange density becomes an 
unreliable approximation to the exact-density as Z decreases 
and aX Z = 1/2 the density is even qualitatively wrong as 
the self-consistent exact-exchange density no longer exists. 
Therefore, we terminate the plot aX Z = 1 where we expect 
the self-consistent exact-exchange density to be an accurate 
representation of the exact density. 



A simple suggestion for f{() which gives both polarized 
and unpolarized limits exactly is 

/(C) = 1 - C (23) 
which is the ( dependence of Uh plus E^. 



IV. DIRACIUM, THE DELTA-FUNCTION 
INTERACTING ANALOG OF HELIUM 

In order to assess the usefulness of this local ap- 
proximation to handle onc-dimcnsional problems, we 
start with perhaps the most difficult test case, a com- 
pletely non-uniform system, diracium. This is the 
one-dimensional analog of Helium with the traditional 
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Coulomb terms replaced by delta-functions. The system 
is described by the Hamiltonian: 



H = 



2dxl 



1 (P 

- Z5{xi) - Z6{X2) + XS{xi - X2) 



and the eigenvalue equation, 



'ixi,X2) 



(24) 



(25) 



whore xi and X2 are the positions of the fermions, a and 
a' are the spin-like labels of the fermions, Z is the magni- 
tude of the external potential, and is an antisymmet- 
ric Fermi wave- function which vanishes as x ^ 00. The 
gromid-statc is a spin-singlct. 

First, we solve the model analytically within the exact- 
exchange approximation. Then, we introduce our local 
density approximation to the cxchangc-corrc^lation en- 
ergy. Finally, we present the exact energy eigenvalues. 

With the exact-exchange functional, Eq. (9), and no 
functional for correlation, the Kohn-sham single-orbital 
equation is 



'2d^x 



<f){x) - Z5{x)(j){x) + X\(j)'^{x)\^(f){x) 



-eKS(l>{x). 

(26) 



Equation (26) can be solved using elementary techniques. 
The resulting eigenvalue is 



eKS = -\[z-\\ 



(27) 



and the Kohn-sham orbital [35] is 
1 



with 



csch 



\x\ 4- a 



a = arccoth 



(28) 



(29) 



This is unbound at Z < ^A. The calculated total energy 
is not just the sum of the KS eigenvalues, rather it is 



E[n] 



e-KS + U[n] +Exc[n] 



j dx VH{x)n{x) — j dx vxc{x)n{x), (30) 



or explicitly, 



v2 ^ 

Ep]xx — — Z -\- — — . 



(31) 



Next, we solve the KS equation using the local density 
approximation to the correlation. The LDA KS equation, 

-\^J>^^) - Z5{x)4>{x) + \\<^\x)\^4>{x) 

+vc{x)<^{x) = -CKsHx) (32) 



with 

vc{x) = 



+ 



2n{x) + d 



n{x) an{x) + b n{x)'^ + dn{x) + e 



(33) 



Equation (32) is solved numerically via a self- consistency 
cycle and the shooting method. 

The exact ground-state energies were obtained previ- 
ously by Rosenthal [36] by transforming to momentum 
space and reducing the problem to the solution of a one- 
dimensional integral equation. While this method con- 
verges quickly to the exact energy eigenvalue, it is not 
well suited to give real space wave-functions and densi- 
ties. Instead, we take the calculated eigenvalue, E, as 
input and reduce the eigenvalue problem, Eq. (24), to a 
differential equation. The differential equation can then 
converted to an integral equation using Green's function 
techniques: 

^{x, y) = - J dx' Ko (y^^/{x - x'y + y^j ^{x', 0) 
+ - j dx' Ko [y^^/x^ + {y- x')^j *(x', 0) 
--J dx' Ko [^/^^/{x - x'Y + {y- x'^j ^{x',x'), 



(34) 



where E is the ground state energy and Kq is the zeroth- 
order modified Bessel-function and is the Green's func- 
tion for this particular two-dimensional equation. By 
setting y equal to x and then 0, we can arrive at a set 
of two coupled integral equations which can be solved 
self-consistently. Note that although the Bessel function 
has a divergence, the integral is finite. Once ^{x, 0) and 
"^{x.x) are known, we can construct the wave-function 
at any point in space. If the inputed E is exact, then 
solution of this equation yields the exact ground-state 
wave- function. 

In Table I, we see that LDA is greatly more accurate 
than the EXX functional. Second order perturbation the- 
ory is only more accurate for Z > 7, but that is at a 
much larger computational cost since the second order 
contribution requires calculation of the entire spectrum 
of unoccupied orbitals. LDA remains bound and gives a 
reasonable result (within a factor of three) even at the 
critical potential strength, Zcrit- In Figure 5, we see that 
the LDA gives a qualitatively correct density for Z = 1 /2 
where the EXX result is no longer even bound. 

There exist standard theorems about the decay of the 
density away from the attractive nuclear potential in the 
three-dimensional Coulomb interacting case [37]. For ex- 
ample, 



n(r) 



-2ar 



(35) 



as r — > 00 where a = V2I. This theorem also holds for 
our one- dimensional model with r = |a;|. There is much 
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interest in critical values of Z at which an atom can no 

longer bind its outermost electron [21 23]. Understanding 
this limit yields information on the existence of negative 
ions. A most interesting question is: As Z ^ Zcxit and 
/ ^ 0, how does the density decay? We study this di- 
rectly in our one-dimensional example by varying A keep- 
ing Z fixed at 1 . For large enough A, the system will ion- 
ize. Note that \crit = ^/Zcrit- Figure 4 shows d\nn/dx 
as A — * 1/Zcrit for Z = 1. 

The second-order perturbation theory result for 
Diracium differs from the high density LDA result. The 
second-order perturbation theory result is: 



Ec = (-3/8 + 2/37r + l/12)\^ 



-o.oygsA-' 



(36) 



The first two terms are the exact contribution to the 
total energy to order A^ [38]. Prom this, we subtract 
the final term which is the exchange contribution via the 
self-consistent density. The high-density LDA correlation 
result is 



Ec = -AVi2 = -0.083333A2 



(37) 



Because the interaction is short-ranged, LDA correctly 
scales to a constant in the high-density limit, and in con- 
trast to three-dimensional Coulomb DFT [39] is highly 
accurate. 

Another interesting quantity to consider as the 
particle-particle interaction grows stronger is the interac- 
tion energy or the expectation value of Vee{xi — Xj). For 
small A, the interaction energy grows in magnitude as A 
grows, but at A w 0.9, this trend reverses. For A > \crit, 
the system is ionized, so that the interaction energy is 
zero. Figure 6 shows that, as A — > Xcrii from below, the 
approach to this discontinuity is linear. This information 
is valuable in studying the approach to ionization, and 
may also be true for real two-electron ions. 



DELTA-FUNCTION INTERACTING 
HOOKE'S ATOM 




X 



FIG. 4: The behavior of the density for Diracium {Z — 1) at 
various interaction strengths, A. We plot dlnn{x)/dx to high- 
light the asymptotic behavior of the density. For A > Xcrit, 
the system is ionized. 
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Another test of this one- dimensional LDA is the analog 



to Hooke's atom [40], 



2dxl 



Id^ 
2dx| 



+ -u> X^ 



^io'^X% + \5{xi 



X2) 



(38) 



where u determines the strength of the harmonic well 
potential. This model has been used to model one- 
dimensional quantum dots [16], and its excitations have 
been studied using time-dependent DFT [41]. The exact 
wave-function is given in terms of Whittaker functions 
and confluent hyper-geometric functions [16]. The total 
energy is 



E = -uj + e 



(39) 



FIG. 5: Comparison of the self-consistent LDA density 
(dashed line) with the exact (solid line) for A = 1 and Z = 0.5. 
This is an extreme case where the exact-exchange no longer 
binds. 



with e obtained from the solution of [16, 41] 

3\ ,„/ e 2' 



2V2wr 



2.^4) 



2u) 



-A. (40) 



In Table II, we see that LDA greatly improves over the 
exact exchange formalism for all values of w. 

In the high-density or weak-coupling (A 0) limit, 
this system behaves similarly to diracium, described 
above. The total energy can be described perturbatively 



E = uj + Aciw^ + A^C2 -I- 



(41) 



7 




FIG. 6: The expectation value of the interaction at various in- 
teraction strengths, A. Beyond A = Acrit (2.6517), the system 
is ionized and the interaction-energy vanishes. 

where ci = 1/V2^ = 0.399 and C2 = (t + /47r = 

—0.110318 with 7 being the Euler constant and ip^°\ the 
zeroth-order polygamma function. The A^-term above is 
the high-density hniit of the correlation energy phis an 
exchange contribution. In DFT, the exchange contribu- 
tion is the first-order contribution in A where A is the 
couphng constant for a fixed fermion density. Here, A 
is the interaction strength, and when it varies so does 
the ground-state density. This means that in a self- 
consistent calculation, exchange will contribute to second 
order in A to the total energy. We saw this in Eq. (31) 
for diracium. The same is true here. Like for diracium, 
LDA is very accurate in the high-density limit. In order 
to properly compared the the LDA correlation energy 
with the exact one, we need to extract the exact density. 

Since Hookc's atom stays bound for arbitrarily weak lo, 
we can see how well our LDA describes the low-density 
limit. The low- density limit is particularly challenging 
for Coulomb-interacting functional [42] because of the 
strong correlation. LDA and generalized gradient ap- 
proximations (GGAs) behave qualitatively correctly but 
often err by as much as a factor of two for the exchange- 
correlation energy. For Hookc's atom, the exact result 
for the total energy in the strong-coupling limit is 

E = 2w. (42) 

The low density limit is related to the large coupling 
constant limit. As A ^ oo, the lowest energy solution 
in the relative coordinate, \xx — X2\, is simply the first 
non-interacting excited state because inserting a node at 
Xi — X2 minimizes the interaction energy. We find that 
LDA is remarkably accurate in this regime, so the errors 
in Coulomb-interacting DFT can be ascribed to the long- 
range interaction. 

The energy in this regime becomes kinetic-like. The 
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TABLE II: Total ground-state energy for the contact- 
interacting Hooke's Atom and errors (in milMHartrees) of var- 
ious approximations with A = 1. Exact is from a numerical 
solution of Eqs. 39 and 40. EXX is exact-exchange. LDA is 
according to the parameterization Eq. 22. 



exact- exchange functional is Hartree-like and will fail to 
capture the proper energetics. The LDA however cancels 
the Hartree-like exchange contributions and is kinetic- 
like. In the low-density regime the density is close to 
uniform locally, so we expect the energy per particle to 
be similar to low-density deltium. This is in fact the case 
and is reflected by the high accuracy of the LDA in the 
low density regime. 



VI. CONCLUSION 

In this paper, we examined a one-dimensional density 
functional theory of contact interacting fermions. We 
noted that exact-exchange is an explicit density func- 
tional and developed a local-density approximation for 
correlation. We have applied these functional success- 
fully to two simple models demonstrating the high accu- 
racy of LDA here. Although LDA is highly accurate in 
these cases, it is not exact. This result is consistent with 
the observation in Ref. [43] that LDA is not exact in the 
short wave-length limit. 

This model interaction and LDA can be used to il- 
lustrate and explore problems in DFT. Examples include 
ground-state symmetry problems in stretched H2 and the 
interacting-fermion one-dimensional solid (a generalized 
Kronig-Penney model) as a model band-gap problem. 
Excited-state Bethe-Ansatz results might be used to de- 
rive a local current-density functional for this interaction. 
This delta-function interaction has already been used in 
scattering problems and in pedagogy, and we hope that 
our local-density correlation functional finds fruitful ap- 
plications in these areas as well. 
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APPENDIX A 

Here, we find the high density hmit of the correlation 
energy per particle for deltium using perturbation theory 
and the diagrammatic approach in momentum space [44] . 
The Fourier transform of the interaction potential is 

V{q) = ^fdx 5{x)e"i^ = j (A-1) 

where L is the length between the periodic boundaries. 
As noted earlier, like spin fermions will not interact via 
the delta-function, this means that only vertices that con- 
nect opposite spins will enter into the perturbative series. 
This is a tremendous simplification as many diagrams will 
cancel. A further simplification is that the interaction is 
independent of the momentum transfer, q. 

To first order, only the Hartree diagram between op- 
posite spins contributes to the total energy. This is the 
first diagram in Figure 7. Evaluation docs not requires 
integration over internal momentum. The loop integrals 
have pre- factors of L^/47r^, and integration over the loops 
results in a factor of 2fci? = wn each, where kp = n/2n 
is the Fermi momentum of non-interacting deltium. The 
symmetry factor of 1/2 is canceled by a sum over the 
two possible pairs of spin. The final energy per particle 
is obtained by dividing by A'', the total particle number. 
We find 

To second order, only one more diagram contributes. 
This is the two-bubble diagram shown second in Figure 
7. From the standard rules of perturbation theory, 

Nec = 



(A-3) 

with \ki + q\> \ and \k2 — g| > 1, fci and ^2 are particle 
momenta, and q is the momentum transfer. The limits of 
integration and constraint inequalities ensure that parti- 
cle have less momentum than the Fermi-momentum, and 



q 




q 



FIG. 7: First and second order contributions to the interac- 
tion energy. Spin labels are omitted but the two loops in each 
diagram must have opposite spins, q is the momentum trans- 
fer, fci and k2 label particle momenta, and ki + q and fe2 — 5 
label hole momenta. 

holes have higher momentum than the Fermi momen- 
tum. Once again, the sum over the two possible spin- 
arrangements cancels the symmetry factor of 1/2. To 
solve Eq. (A-3) exactly, we rescale as follows: q = kpx, 
ki = kpy, and k2 = kpz. After some algebra, we find 
the correlation energy per particle: 




(A-4) 



using the quadrature result, 

1 = 2 1 dx dy dz —, r- 

A J-\ J-i x{x + y-z) 

1-2 j-l r-l+x ^ 27r2 

+2 dx dy dz — — = — . 

Jo Ji-x J-1 x{x + y-z) 3 

(A-5) 
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